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We study numerically the dependence of heat transport on the maximum velocity and shear rate 
of physical circulating flows, which are prescribed to have the key characteristics of the large-scale 
mean flow observed in turbulent convection. When the side-boundary thermal layer is thinner than 
the viscous boundary layer, the Nusselt number (Nu), which measures the heat transport, scales 
with the normalized shear rate to an exponent 1/3. On the other hand, when the side-boundary 
thermal layer is thicker, the dependence of Nu on the Peclet number, which measures the maximum 
velocity, or the normalized shear rate when the viscous boundary layer thickness is fixed, is generally 
not a power law. Scaling behavior is obtained only in an asymptotic regime. The relevance of our 
results to the problem of heat transport in turbulent convection is also discussed. 

PACS numbers 47.27.-1,44.20. +b 



I. INTRODUCTION 

Turbulent Rayleigh-Benard convection has been a sys- 
tem of much research interest. The system consists of a 
closed cell of fluid which is heated from below and cooled 
from above. When the applied temperature difference is 
large enough, the fluid moves and convection occurs. The 
flow state is characterized by the geometry of the cell and 
two dimensionless control parameters: the Rayleigh num- 
ber Ra, which measures how much the fluid is driven and 
the Prandtl number Pr, which is the ratio of the diffusiv- 
ities of momentum and heat of the fluid. The two param- 
eters are defined by Ra = agAL^ /{vk), and Pr = v/k, 
where A is the maintained temperature difference, L is 
the height of the cell, g the acceleration due to gravity, 
and a, v, and k are respectively the volume expansion 
coefficient, kinematic viscosity and thermal diffusivity of 
the fluid. When Ra is sufficiently large, the convection 
becomes turbulent. 

Besides the issue of the statistical characteristics of the 
velocity and temperature fluctuations, it is of interest to 
understand the heat transport by the fluid, which is an 
overall response of the system. The heat transport is usu- 
ally expressed as the dimensionless Nusselt number Nu, 
which is the ratio of the measured heat flux to the heat 
transported were there only conduction. Before the on- 
set of convection, heat is transported only by conduction 
and Nu is identically equal to one. When convection oc- 
curs, heat is more effectively transported by the fluid due 
to its motion and Nu increases from 1. A major question 
is then to understand how Nu depends on Ra and Pr. 

The work of Libchaber and coworkers on turbulent 
convection in low temperature helium gas ||l|,^ showed 
that Nu has a simple power-law dependence on Ra: 



Nu - Ra 







(1) 



and the exponent /3 is almost equal to 2/7, which is dif- 



ferent from 1/3, the value that marginal stability argu- 
ments would give. This result led to the development 
of several theories |2|,^-^ which all give (3 = 2/7 but 
were based on rather different physical assumptions. In 
particular, the Chicago mixing-zone model Q| empha- 
sizes the heat transport by the thermal plumes, the co- 
herent structures observed in turbulent convection while 
the theory by Shraiman and Siggia focuses on the 
effect of the shear of the large-scale mean flow on the 
heat transport (see e.g. Ref. for a review). Recent 
experimental results appeared to further complicate the 
situation. Niemela et al. ||^ reported a value of (3 close 
to 0.31 for measurements in low temperature helium gas 
that cover a much larger range of Ra, from 10^ to 10^^. 
Xu et al. |9| studied turbulent convection in acetone in 
several experimental cells of different aspect ratios and 
concluded that there is no signiflcant range of Ra over 
which the scaling behavior Eq. (|^) holds. Furthermore, 
they showed that for Ra > 10^, the dependence of Nu 
on Ra can be represented by a combination of two power 
laws, which is consistent with a recent theory by Gross- 
mann and Lohse ||l^,|l^. In Grossmann and Lohse's the- 
ory, the viscous and thermal dissipation were decomposed 
into their bulk and boundary-layer contributions, and ten 
asymptotic regimes were obtained |ll| . 

Another interesting feature observed in turbulent con- 
vection is the presence of a persistent large-scale mean 
flow which spans the whole experimental cell ||l^. The 
maximum mean velocity of the flow was also found to 
scale as Ra to about 1/2 [|l3|. The presence of a large- 
scale flow naturally induces an interaction between the 
top and bottom thermal boundary layers. Such an in- 
teraction was taken to be absent in the marginal stabil- 
ity arguments. One obvious effect of the velocity held, 
which satisfies the no-slip boundary condition, is that it 
produces a shear near the boundaries, which was first 
studied in Ref. Istl. 
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In convection, the equations of motion are: 

u ■ Vw = — Vp + vV^u + gaTz 
u-VT = kV^T 
V • M = 



du 

'm 



dT 
~dt 



(2a) 
(2b) 
(2c) 



where u is the velocity field, p the pressure divided by 
density and T the temperature field, while z is the unit 
vector in the vertical direction. The velocity and tem- 
perature are thus coupled dynamically in a complicated 
fashion and have to be solved together. Physically, the 
velocity field is driven by the applied temperature dif- 
ference. The flow in turn determines the temperature 
profile, and thus the heat transport, in a self-consistent 
manner. 

To gain insights of the problem of heat transport in tur- 
bulent convection, we have turned to the simpler problem 
of heat transport by prescribed velocity fields that have 
features of the large-scale mean flow observed. The large- 
scale flow has two dominant features: (i) it is a circulat- 
ing flow that spans the whole experimental cell and (ii) it 
generates a shear near the boundaries. In an earlier pa- 
per, Ching and Lo studied separately these two features 
and their effects on the heat transport They found 
that Nu scales with the Peclet number that measures the 
maximum velocity to an exponent 1/2 for a purely circu- 
lating flow and scales with the normalized shear rate to 
an exponent 1/3 for a pure shear. 

Pure circulating or pure shear flows are, however, not 
physical fluid flows within a closed box. In this paper, we 
continue along these lines of thought and study the heat 
transport by physical velocity fields in a unit square cell 
that have both the circulating and the shear-generating 
features. We focus on the dependence of the heat trans- 
port on the maximum velocity and the shear rate of the 
flows. We first formulate our problem in Sec. II. Then we 
present our numerical results and discuss how these re- 
sults can be understood in Sec. III. In Sec. IV, we further 
discuss how these results are relevant to the understand- 
ing of heat transport in turbulent convection. Finally, we 
end the paper with a summary and conclusions in Sec. 
V. 



II. THE PROBLEM 

We solve the steady-state advection-diffusion equation 

u{x,y)-VT{x,y) = >iV^T{x,y) (3) 

for a prescribed incompressible velocity field u{x, y) in a 
unit square cell: < a; < L and < y < L. A tempera- 
ture difference of A is applied across the y-direction while 
no heat conduction is allowed across the x-direction. 
That is, the temperature field T{x,y) satisfies the fol- 
lowing boundary conditions: 



T(x,y-0) = A; 

dT 

-(x = 0,y) 



T{x,y = L) = 
dT 

{x = L,y) 



dx ' dx 

For the velocity field, we take 

Ux{x,y) = f{x)f'{y) 
Uy{x,y) = ~f'{x)f{y) 



(4) 
(5) 



(6) 
(7) 



where / is some function oix = x/Lovy = y/L and ' 
is its derivative with respect to the argument. Thus, u is 
incompressible and separable in x and y. To satisfy the 
no-slip boundary condition, we require 

/(O) = /(I) = no) = /'(I) = (8) 

Moreover, for u to be a circulating flow, we require 

/(5) = /(l-i) 0<i<l (9) 

such that Ux and Uy are antisymmetric about y — L/2 
and X — L/2 respectively. We have studied two forms 
of /. The first form is algebraic. When / is algebraic, 
Eqs. (|) and (|) imply that /(i) = i'^{i- lfh{i) with 
h{x) — h{\ ~ x). Thus we choose 



f{x) = x^i - lY{ai + b)[a{l -S:)+b] 



(10) 



where a and b are positive constants. Equation ([lo| ) is the 
simplest algebraic form that allows us to change both the 
circulating strength and the shear rate of the flow (see 
below). The other is an exponential form of /: 

f{i) = c(l - e-'^*)2[l - e-'^^i-*)]^ (11) 

where c and k are positive constants. For fc <gC 1, the ex- 
ponential form reduces to the algebraic form with a = 
and b^ — ck^ . When k is large, the velocity decays ex- 
ponentially towards the center of the cell. We show the 
two velocity fields in Fig. |l|. 
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FIG. 1. (a) The velocity field with an alegbraic / given 
by Eq. ( p^ with a — b = 1. The size of the arrow indicates 
the relative magnitude of the velocity. The contour C that 
encloses the lower half of the unit square cell is also shown, 
(b) The velocity field with an exponential / given by Eq. ( |ll| ) 
with c = 1 and k = Q. 
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We characterize the velocity field by its maximum ve- 
locity Mo and shear rate 7, which are defined by 



1.0 



Mo = max Uxix = —,y) 

0<y<i 2 



max -Uy{x,y= -) (12) 
o<x<4^ 2, 



dy 



L 



V) 



y=o 



du. 
dx 



(13) 



x=0 



By varying the parameters a and 5 or c and fc, we can 
vary mq and 7 of the velocity fields. 

From the solved temperature field, we can calculate 
the heat transport, which is the heat conducted across 
the boundary y = 0. Nu is therefore the ratio of the aver- 
age of the magnitude of the vertical temperature gradient 
over the boundary y = to A/L: 



Nu 



dT 
dy 



A 
1 



(14) 
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FIG. 2. Typical vertical and horizontal temperature pro- 
files T{x = L/2,y)/A (solid) and T{L - x,y = L/2)/A (dot- 
ted) as a function of y/L and x/L respectively. 
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Here (. . .) is the average over x from to L. Our goal is 
to study the dependence of Nu on the Peclct number Pe 
= uqL/k and the normalized shear rate 7 = jL'^/k. 



III. RESULTS AND DISCUSSIONS 
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We numerically solve T{x,y) for the two forms of /. 
In Fig. ^, we show the vertical and horizontal temper- 
ature profiles T{x = L/2,y) and T{L - x,y = L/2). 
As was reported in turbulent convection experiments, 
the applied temperature difference concentrates in two 
narrow regions near the 'bottom' and 'top' boundaries 
y = Q and y = L respectively. Interestingly, the cir- 
culation also induces a temperature difference between 
the two 'side' boundaries a; = and x = L. For flows 
that are anticlockwise (see Fig. the side x = is 
colder than the side x = L. The horizontal tempera- 
ture profile T{L — x,y — L/2) resembles the vertical one 
T{x — L/2,y) in that it is almost constant = A/2 ex- 
cept for two small regions near the sides. In these two 
small regions, it is approximately linear in x except at 
the boundaries where the horizontal gradient vanishes. 
Hence, we approximate T(x, y = L/2) by 



T{x,y^L/2) 
A/2 




a 1 



ail 



L—x ^ 



< a; < ^ 
l<x<L-i 
L-e<x< L 



(15) 



where a is about 0.5 as shown in Fig. ^ 
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FIG. 3. The dependence of a on 7 for various values 
of A„: \y/L = 0.0167 (circles), A„/L = 0.021 (trian- 
gles), \y/L = 0.028 (stars), \y/L = 0.070 (squares), 
\v/L = 0.088 (plusses), X^/L = 0.096 (diamonds), 
\v/L = 0.108 (crosses), and X^/L = 0.130 (inverted trian- 
gles). 

Thus, there are also two thermal boundary layers of 
thickness £ at the 'side' boundaries. Furthermore, we 
find that Nu is given by L/£ up to a factor d: 



Nu 



.L 



(16) 



where d is weakly independent of 7 and approaches 0.8 as 
the viscous boundary layer thickness A„ = uq /j increases 
to 0.13 L as shown in Fig. |. Defining At = L/(2Nu) 
to be the average thickness of the thermal boundary 
layer at the top and bottom boundaries, we see that the 
£ = 2d\T- Thus, the side- wall thermal boundary layers 
are thicker than the top and bottom thermal boundary 
layers. 
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FIG. 4. The dependence of d on 7 for various values of 
A„ (same symbols as in Fig. 3). 



We find that the functional form of Nu(Pe,7) depends 
crucially on the relative sizes of At, and £. For £ < X^, 
Nu scales vifith the shear rate: 



Nu = Aj 



^1/3 



for f < A„ 



(17) 



and the coefficient A approaches 0.36 as A„ increases to 
0.13 L as shown in Fig. ^. For A„ < £, we see in Fig. ^ 
that the dependence of Nu on Pe or 7 for fixed Xy is not 
a power law. For a short range of Pe or 7, the depen- 
dence might be described by an effective power law but 
the value of the effective exponent would depend on the 
range of Pe or 7 fitted. 

We shall understand these results in the following. 
Since the velocity field is incompressible, Eq. (||) implies 



uT - KVTj ■hdl = 



(18) 



for any closed curve C in the two-dimensional domain 
where h is an outward normal. Since we find that \dT / dy\ 
almost vanishes along y = i/2, we choose C to enclose 
the lower half of the unit cell (see Fig. |l|a). Together 
with Eq. (^ and the no-slip boundary condition, Eq. 
( p^ ) then implies that 

^n^^{uy{x,v = ^)T{x,v = ^)) (19) 

Thus, we can estimate Nu by the heat transported across 
y = L/2. Using the antisymmetry of Uy about x — L/2 
and Eq. (15), we get 



{uy{x,y^ ^)T{x,y ■ 
- / Uy{x,y= -) 



L 



)) 



' Jo 
aA 



L, A' 

nx,y^-)-- 



dx 



(1 - -)uy{x,y = —)dx 



(20) 



In the derivation of Eq. (|2C|), we have assumed that 
the parameters of the velocity field are chosen such that 
£ < -L/2 as in physical situations. The heat transported 
across y — L/2 is, therefore, contributed mainly by two 
'jets' of colder and hotter fluids moving down and up 
respectively along the two boundaries x — Q and x — L. 



100 



10 



0.375 








□ 


□ □□ □ 




^ 0.365 
0.355 


v"" VW 


i+ +++ 




10= 


10° 




y 




+ + 










V 









10" 



10° 



FIG. 5. The dependence of Nu on 7 when I < \v for var- 
ious values of A„ (same symbols as in Fig. In the inset, 
Nu7~^^^ is plotted versus 7. 
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FIG. 6. The dependence of Nu on 7 when L/2 > ^ > A„ 
for various values of A„ (same symbols as in Fig. We also 
compare our numerical results to Eq. ( p^ with X^jL — 1/60, 
a = 0.6 and d = 0.9 (solid line), and good agreement is found 
for i < O.IL (data points on the right of the dashed line). 



Using Eq. (|T^ , we approximate Uy by a linear function 
in x for a; < Ai,: 



Uy{x, — ) — —JX for X < Xy 



(21) 



Hence, Eqs. (|l|) and give 
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Nu « f ^ j 7^/3 for £ < A„ (22) 

In Fig. 1^, we compare our estimated values of the co- 
efficient with the computed values of A. It 
can be seen that the two values agree within an error of 
7%. Moreover, the agreement is better for larger Xy, as 
expected since the linear approximation Eq. IGV) works 
better ior < x < £ for larger A^. 



depend on the range of Pe or 7 fitted. In Fig. we com- 
pare the numerical results to Eq. (^ ) with \y/L = 1/60, 
a = 0.6 and d — 0.9. The agreement is not too bad for 
Xy < £ < O.IL given that our approximation of Uy being 
constant beyond the shear layer is rather crude. 



IV. RELEVANCE TO TURBULENT 
CONVECTION 
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FIG. 7. Comparison of the estimated values of the co- 
efficient {ad^/Gy^^ (circles) with the computed values of 
A (squares). 



For £ > Xy, we need Uy beyond the region where a lin- 
ear approximation holds. As a first approximation, we 
take 



7a; 

. Uf) 



X < Xy 

Xv < X < 



(23) 



That is, we approximate the large-scale flow by a shear 
near the boundaries then followed by a band of circula- 
tion at the maximum velocity. As Uy has to decay to zero 
towards the center of the cell, we expect the approxima- 
tion Eq. to work only for Xy < £ ^ L/2. 

Using Eqs. ([l9|), (|o|) and (^, we get a quadratic 
equation for Nu. Solving which gives 



Nu = 



2adPe - (aVl2)Pe''7-2 - (a/2)Pe^7 



(24) 



for L/2 ^ £ > Xy. Thus, Nu does not have a power- 
law dependence on Pe nor 7. Physically, this is because 
we generally cannot neglect the effect of the shear even 
when the edge of the side-wall thermal boundary layer 
is located at the band of maximum velocity of the large- 
scale flow. It is only in the limit Pe'^/^/7 <^ 1 that Nu 
~ Pgi/2_ Away from this asymptotic regime, Nu might 
be represented by an effective power law on Pe or 7 for 
fixed Xy but the value of the effective exponent would 



In turbulent convection, the heat transport is due both 
to the large-scale mean flow and the fluctuating part of 
the velocity field. Our present work provides insights 
only to the heat transport by the large-scale mean flow. 
It is illuminating to see what results would be inferred if 
we neglect the effect of the fluctuating part of the velocity 
field. 

Depending on the type of the viscous boundary layer, 
the strength and the shear rate would be related to each 
other. For example, if we follow Grossmann and Lohse 
to assume that viscou s boundary layer is of Blasius type 
| p5[ : Xy/ L ^ -^/Pc/Pr for moderate values of Pr jl^ and 
that Xy ^ L for very large values of Pr , then 
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Pe^/^Pr"^/^ moderate Pr 
Pe very large Pr 



(25) 



Our results Eqs. (|2^) and (^4|) would thus give Nu as a 
function of Pe and Pr: 



Nu 



1/3 



^r^' Pel/3 



Pe^/^Pr-i/^ moderate Pr 
very large Pr 
for £ < Xy and 

Nu « g(Pr)Pei/^ 
for L/2 > £ > Xy with 

y^ad- a^Pr/U ~ aV?'^ /2 



g(Pr) 



2[1 - aPr/(6d)] 



(26) 



(27) 



(28) 



We emphasize that the non-power-law dependence on Pr 
for £ > Xy echoes that the effect of the shear cannot 
generally be neglected even when the edge of the side- 
wall thermal boundary layer is located at the band of 
the maximum velocity of the large-scale flow. 

Next, we make use of a rigorous relation between the 
viscous dissipation and Nu and Ra, which is derivable 
from the equations of motion ||^,|^ : 



(([9.u,(x,t)]2)) = — (Nu-l)Ra 



(29) 



where ((. . .)) is an average over space x and time t. Ne- 
glecting the contribution from the fluctuating part of the 
velocity field to the viscous dissipation, we get 



7Pe NuRa 



(30) 
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for Nu > 1. Using Eqs. (g5|), (g6|), (g7|) and (|3C|), we 
finally get 



Nu ^ Pr-l/l^Ral/^ Pe ^ Pri/'^Rai/^ 



Nu ~ Ra 



1/5 



Pe ~ Ra' 



3/5 



(31) 



(32) 



for £ < Xy respectively for moderate and very large values 
of Pr, and 



Nu - 5(Pr)5/4pri/8Rai/4 
Pe ^ g(Pr)i/2pri/4Rai/2 



(33) 



for L/2 > ^ > A„. The results Eqs. (|3l|), (g|), and 
(^), except for the non-power-law dependence on Pr 
in Eq. (p3|), resemble those obtained by Grossmann and 
Lohse respectively for regimes //, and /„, in 
which they took the boundary layer to dominate both 
the viscous and thermal dissipation. 



not a power law and scaling behavior is obtained only in 
the asymptotic regime Pe'^/^/7 <C 1. 

We have further discussed how our results are relevant 
to the problem of heat transport in turbulent convection. 
In turbulent convection, the heat transport is due both 
to the large-scale mean flow and the fluctuating part of 
the velocity field. Neglecting the effect of the fluctuating 
part of the velocity field, our results lead to results re- 
sembling those obtained by Grossmann and Lohse [pO 11 
when they took the boundary layer to dominate both the 
viscous and thermal dissipation. It is not surprising that 
the boundary layer dominating the viscous dissipation is 
the same as the large-scale mean flow dominating the vis- 
cous dissipation since the gradient of the large-scale mean 
flow, which contributes to the viscous dissipation, con- 
centrates in the boundary. Our finding thus suggests that 
whether the boundary layer or the bulk dominates the 
thermal dissipation is physically equivalent to whether 
the large-scale mean flow or the fluctuating part of the 
velocity field dominates the heat transport. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the heat transport by 
physical fluid flows whose velocity fields are prescribed to 
have both the two key characteristics of the large-scale 
mean flow observed in turbulent convection. The velocity 
fields that we have chosen are separable and incompress- 
ible circulating fiows in a unit square cell. Satisfying the 
no-slip boundary condition, they also generate a shear 
near the boundaries. Overall, they are approximately a 
shear near the boundaries, almost constant with the max- 
imum strength of circulation for a finite band, and then 
a decay towards the center of the cell. We focus on the 
functional dependence of Nu on Pe measuring the max- 
imums strength of circulation and the normalized shear 
rate 7 that characterize the velocity fields. 

We have shown that Nu can be estimated by the heat 
transported across the mid-height of the unit square cell, 
which is in turn contributed mainly by two jets of hot- 
ter and colder fluids moving up and down the two sides. 
These two jets are confined to two narrow regions. That 
is, the velocity field also induces two thermal boundary 
layers at the side boundaries. These side-boundary ther- 
mal layers are thicker than those at the top and bottom 
boundaries. It is then clear that the functional form of 
Nu(Pe,7) depends crucially on the relative sizes of the 
viscous boundary layer thickness and the thickness 
of the thermal boundary layers at the side boundaries i. 
When ^ < A^, that is, the edge of the side- wall thermal 
boundary layer falls within the shear region of the large- 
scale flow, Nu scales with 7 to 1/3 and is very weakly 
dependent on A„ or Pe. When £ > X^ and the edge 
of the side-wall thermal boundary layer falls within the 
band of circulation with maximum strength, there is still 
contribution to the heat transport by the shear and Nu 
depends on both Pe and 7. The dependence is generally 
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